Fragmentation of spin-dipole strength in Zr and Pb 
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Abstract 

An extension of time-dependent covariant density functional theory that in- 
cludes particle-vibration coupling is applied to the charge-exchange chan- 
nel. Spin-dipole excitation spectra are calculated an compared to available 
data for ^°Zr and ^"^Pb. A significant fragmentation is found for all three 
angular-momentum components of the spin-dipole strength as a result of 
particle-vibration coupling, as well as a shift of a portion of the strength to 
higher energy. A high-energy tail is formed in the strength distribution that 
linearly decreases with energy. Using a model-independent sum rule, the cor- 
responding neutron skin thickness is estimated and shown to be consistent 
with values obtained at the mean-field level. 

Keywords: covariant density functional theory, particle-vibration coupling, 
nuclear charge-exchange excitations 

Spin-isospin excitations present a very active research topic both in nu- 
clear structure and nuclear astrophysics. In particular, a detailed knowledge 
of the Gamow- Teller resonance, a collective oscillation of excess neutrons that 
coherently change the direction of their spin and isospin without changing 
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their orbital motion, is essential for understanding weak nuclear reactions in- 
volved in the process of nucleosynthesis, i.e. /3-decay, electron and neutrino 
capture. Moreover, it has been shown that spin-dipole charge-exchange ex- 
citations, made up of three components with angular momentum and parity 
J'^ = 0^, 1~ and 2", can significantly contribute to the total reaction rates 
and even compete with the contribution of Gamow- Teller transitions l|,l2|, |3|. 

The spin-dipole strength can also provide information on basic proper- 
ties of finite nuclei. The thickness of the neutron skin has been shown to 



constrain the neutron equation of state 



^, and is also correlated with the 



nuclear symmetry energy . While accurate data on the charge distribution 
in nuclei have been obtained by elastic electron scattering, our knowledge of 
neutron distribution comes primarily from hadron scattering, and the results 
are markedly model-dependent. Indirect methods for determining the neu- 
tron skin thickness have been proposed, based on energy differences between 
the Gamow- Teller and isobaric analogue resonances^, and using the model- 
independent sum rule for the spin-dipole resonance [7] . Two recent (p^n) and 
{n,p) measurements of the spin-dipole response of ^"Zr and ^osp]-, g| |g|j j^^ve 
also prompted new theoretical studies, in particular investigations based on 
the random phase approximation [ll|, Q, 13|. 

Both measurements ^, ol show a high-energy tail in the spin-dipole 
strength distribution that cannot be described by the simple one-particle 
- one-hole (Iplh) random phase approximation (RPA). Previous attempts 
to extend this framework using the 2p2h RPA were based on a non self- 
consistent approach that employs a phenomenological Woods-Saxon poten- 



tial to obtain the ground-state wave functions 



14| . In this Letter we in- 
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troduce a charge-exchange version of the particle-vibration couphng model 
based on time-dependent covariant density functional theory, and apply it to 
an analysis of spin-dipole strength distributions in ^°Zr and ^'^^Pb. 

The basic quantity that describes small-amplitude motion of an even-even 
nucleus in an external field with frequency u is its response function R{uj) 
[isl - It is obtained as a solution of the Bethe-Salpeter equation: 

R{uj) = R\uj) + lf{uj)W{uj)R{uj), (1) 

where R^{uj) is the propagator of two uncorrelated quasiparticles in the static 
mean field, and the second term includes the in-medium nucleon-nucleon 
interaction W{u). The two-body interaction W{u) contains static terms 
and a frequency-dependent term: 

W{uj) = Vp + V^ + Vt^ + - $(0). (2) 

Vp and VJr represent the finite-range p-meson and vr-meson exchange inter- 
actions, respectively. They are derived from the covariant energy density 
functional and read: 

1^(1,2) = (72^ir2(/37'')i(/37M)2^p(ri,r2), 

V;(l,2) = (^)Vir2(SiVi)(5]2V2)/^.(ri,r2), (3) 

where gp and /^r are the coupling strengths, Dg a nd D^, are the meson propa- 



gators and S is the generalized Pauli matrix IGj. The derivative type of the 



pion-nucleon coupling necessitates the inclusion of the zero-range Landau- 
Migdal terni Vl'^^ , which accounts for the contact part of the nucleon-nucleon 
interaction jl6|. $(0;) describes the coupling of the quasiparticles to vibra- 
tions (phonons) generated by coherent nucleonic motion. In the quasiparticle 



time blocking approximation (QTBA) [17] it can be written in the following 
operator form: 

H^) = ^^"^'^^ - V^m)9^r:i\ (4) 

where the index m enumerates vibrational modes with frequencies Um and 
coupling amplitude matrices gm\ and the index rj = ±1 denotes forward and 
backward components. In Eq. ([2]) the term $(0) is subtracted to remove the 
effect of double counting the phonon coupling, because the parameters of the 
density functional have been adjusted to ground-state data and, therefore, 
already include essential static phonon contributions. The energy- dependent 

^mentation of nuclear spectra and 



18| 



effective interaction of Eq. @ leads to fra 
determines the width of giant resonances 
The strength function S{u): 

S{E, A) = --Im{P^R{E + iA)P), (5) 

TT 

yields the spectral distribution of the nuclear response in a given external 
field P. The field operators for charge-exchange spin-dipole transitions read: 

^i = E^w['^w®^i«]A^±w, (6) 

i 

where t± denotes the isospin raising and lowering operators. 

The phonon coupling terms augment the RPA spectrum with additional 
p-h®phonon components that generally lead to significant fragmentation of 
giant resonances 18| . In Fig. [T] we compare the spin-dipole strength distri- 
bution in ^°Zr calculated using the relativistic random phase approximation 
(RRPA), and the relativistic time-blocking approximation (RTBA). In both 
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models the NL3 [19| relativistic mean-field effective interaction has consis- 
tently been used for the calculation of the mean-field ground-state, the RPA 
phonons, and the spin-dipole charge-exchange excitations. 

The prominent RRPA peaks (dashed curves) disappear when particle- 
vibration coupling is included in the RTBA (solid curves). In the t„ channel, 
for instance, only a broad resonance remains with the peak at 23.5 MeV, in 
very good agreement with available data 8|]. The three angular- momentum 
components do not, however, follow the energy hierarchy E(2^) < E{1^) < 



E{0~) predicted by recent Skyrme-RPA calculations 10|, While the 

RTBA predicts the 2~ component to be the lowest one with the centroid 
energy ml / itlq = 25.4 MeV, 0~ is found to be lower than the 1~ component, 
with centroid energies at 29.4 MeV and 32.6 MeV, respectively. This is 
probably due to the fact that the exchange terms are neglected in the mean- 
field calculation. Namely, as shown in a recent study |12)], a fully consistent 
relativistic Hartree-Fock (RHF) + RPA calculation yields E{2') < E{1') < 
E{0^) for the excitation energies of spin-dipole components. 

The inclusion of particle- vibration coupling leads to a shift of the strength 
to higher excitation energies. A high-energy tail is formed in the region above 
30 MeV where the strength decreases almost linearly with increasing energy, 
in close agreement with experimental results. In contrast, the RRPA strength 
decreases more rapidly above 30 MeV, and becomes 5 to 10 times smaller 
than the experimental strength above 40 MeV. 

In the t+ channel the two dominant peaks predicted by the RRPA merge 
into a single broad structure that extends up to approximately 15 MeV ex- 
citation energy. The tail at higher energies decreases approximately linearly 
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Figure 1: Spin-dipole strength distributions for the t_ (upper panel) and the t+ (lower 
panel) channels in ^'^Zr. On the horizontal axis the excitation energy is plotted with 
respect to the ground state of ^''Zr. The solid black curve represents the sum of the 
strength distributions of the 0^, 1^ and 2^ components, calculated in the RTBA that 
includes particle- vibration coupling. The dashed black curves denotes the corresponding 
RRPA strength function. In both cases the imaginary part of the energy is set to A = 1 
MeV. The experimental results denoted by full circles are from Ref. jsj. 

with increasing energy. One might notice a very good agreement with data, 
except in the low-energy region below 5 MeV, where both the RRPA and the 
RTBA predict spin-dipole strength, originating predominantly from the 2~ 
component, that is considerably larger than the measured distribution. 

The results for the spin-dipole strength in ^°'^Pb are shown in Fig. |2j 
The data exhibit a broad asymmetric resonance centered at 25 MeV, and an 
additional small peak at approximately 6 MeV. The RTBA results reproduce 
these structures, even though the calculated width of the main resonance is 
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slightly larger than the empirical value. As in the previous case, a portion of 
the strength is shifted to higher energies by the inclusion of particle- vibration 
coupling, in very good agreement with data above 35 MeV. Obviously in this 
region the RRPA strength distribution decreases much faster with energy 
compared to RTBA. The ordering of the angular-momentum components of 
the strength is the same as in the case of ^°Zr: E{2^) < E{0^) < E{1^). 
The lower panel of Fig. [2] displays the distributions of the SD+ strength. In 
this case the strength is concentrated in a single peak centered at 7.5 MeV. 
Relatively little fragmentation is obtained in comparison with the RRPA 
results, even though some strength is shifted to higher excitation energy in 
the RTBA. These findings are supported by the available data from Ref. 20|, 



where a single peak has been observed at approximately 7.5 MeV excitation 
energy. 

In the case of ^°'^Pb data are also available for each component of the spin- 
dipole strength [oj. In Fig. Elwe display a comparison between the RRPA, 
RTBA and the experimental results for J'^ = 0~, 1~ and 2^. The RRPA 
predicts that the strength of the 0~ component is concentrated in a single 
peak at the excitation energy of ~ 28 MeV with respect to the ground state 
of ^°^Pb. Particle-vibration coupling induces fragmentation and spreading 
of this strength, but the basic structure of the distribution is not altered. 
Obviously this does not completely agree with the experimental results. For 
the 1~ components the main peak is centered around 30 MeV, whereas the 
experiment places it around 23 MeV. The opposite situation occurs for the 2~ 
component, for which the calculated distribution is in qualitative agreement 
with experiment, even though the centroid of the main peak is calculated 
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Figure 2: Same as described in the caption to Fig. [T]but for ^"*Pb. The data are from 

Rcf. lal. 



few MeV below the measured resonance. In Ref. 13| the results have been 
brought in agreement with experiment by the inclusion of tensor correla- 
tions in the Skyrme energy density functional. These correlations exhibit 
a mult ipole- dependent effect on spin-dipole excitations. It is interesting to 
note that in this case experimental results indicate that the 1~ component 
of the spin-dipole resonance is actually below the 2~ component. 

Using the calculated strength distributions, we have also determined the 
sum-rule values for the spin-dipole response of ^^Zt and ^"^^Pb. The model 
independent sum rule relates the spin-dipole strength to the neutron and 



proton ground-state radii 



211: 



S'--Sl^'-^{N{r\-Z{r')_;), (7) 
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Figure 3: Spin-dipole strength distributions in ^°^Pb for the J'^ = 0~,1~ and 2~ com- 
ponents. The black dashed curves are the RRPA results, and the solid red curves denote 
results obtained with the particle- vibration coupling model. Data are from Ref. 
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Table 1: Spin-dipole sum-rule values, proton and neutron rms radii, and thickness of 
the neutron skin of ^'^Zi and ^°*Pb. The proton radii correspond to the self-consistent 
ground-state distributions, whereas the neutron radii are calculated using the sum rule 
of Eq. ([7]). Sum rules are given in units of fm^, and radii in fm. Error estimates in 
the experimental value of the sum rule arise from statistical, systematic and multipole 
decomposition uncertainties, respectively. 





902r 


208p]jj 


S- - S+ (g.s.) 


160.925 


1222.044 


S- — 5+ (calc.) 


160.963 


1213.562 


S- — S+ (exp.) 


148 ± 6 ± 7 ± 7 [8J 




{r^)p (g.s.) 


4.193 


5.459 


\/(^^)n (calc.) 


4.308 


5.731 


Vi^^)n (exp.) 


4.26 ±0.04 [8j 




6np (calc.) 


0.115 


0.272 


6np (exp.) 


0.07 ± 0.04 [8] 


0.083 < 5np < 0.111 [25] 
0.156^°:°^? [22 
0.19 ± 0.09 [23] 



where denotes the total SD strength in the t± channel for angular mo- 
mentum A. S± will denote the sum of spin-dipole strengths of the three 
components. The neutron skin thickness Snp is defined as the difference of 
neutron and proton rms radii: 

Snp=^/^r.-^p■ (8) 

The results obtained in the present study are summarized in Table [1] The 
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first and second rows of tlie table give tlie values of the sum rule obtained 
employing the self-consistent ground-state mean-field solutions and the cal- 
culated SD strengths, respectively. The experimental value of the sum rule 
for ^°Zr jsl is shown in the third row. The difference between the two the- 
oretical results is very small because the model is self-consistent, but they 
both overestimate the experimental value. The radii of proton distributions 
in the fourth row are extracted from the self-consistent ground-state densi- 
ties, and these values are in excellent agreement with data 22|]. In the fifth 
and sixth rows we include the rms radii of neutron distributions in ^°Zr and 
^°^Pb, calculated from Eq. ([7]) using the proton ground-state radii and the 
calculated SD strength distributions, and the experimental value for ^°Zr, 
respectively. Finally, the calculated and experimental values for the neutron 
skin thickness are given in the last two rows of Table [H 

For ^''Zr the calculated sum rule is 9% larger than the measured value. 



This leads to a neutron rms radius ^/Jr^)^ = 4.308 fm, and neutron skin 
thickness 6np = 0.115 fm, both at the upper limit of the experimental error 
bars. We note that the relativistic Hartree-Fock -|- RPA calculation predicts 



the neutron skin thickness 5„p = 0.092 fm 



12 [, while Skyrme-based results 



io|. 



range from 6np = 0.055 fm to 6np = 0.106 fm 

There are no experimental values of the total spin-dipole strength in 
^°^Pb, but several measurements of the neutron skin thickness have been 
reported. A comparison of the measured cross section for the isoscalar gi- 
ant dipole resonance and the DWBA calculation yielded the neutron skin 

23|. Microscopic optical potential analyses 



thickness 6np = 0.19 ± 0.09 fm 



of intermediate energy elastic proton scattering give 5„p 0.17 fm 



2^, and 



11 



0.083 fm < 5np < 0.111 fm [25!]. Using the correlation between neutron skin 
thickness and the isovector dipole polarizabihty obtained with the Skyrme 



functional 



26| . the very recent results on polarized proton inelastic scatter- 



ing at forward angles yield the value 6np = 0.1561o;q2i fm [27]. The value 
obtained from the calculated SD sum rule in the present study is 5„p = 0.272 
fn,, —ably la.ge. than the en,pi.ica, value. The RHF . RPA .odel 
predicts 6np = 0.234 fm |12|, and values obtained from various Skyrme-based 
models range between 6np = 0.125 fm and Snp = 0.228 fm |10]. The relatively 
large neutron skin thickness that we have obtained in the present study is 



19|, characterized by a 



peculiar to the relativistic effective interaction NL3 
large asymmetry energy. Employing one of the modern relativistic function- 
al with non-linear effective interactions in the isovector channel and lower 
asymmetry energy as, for instance, DD-PCl [28], the thickness of ground- 
state neutron skin is calculated: 6np = 0.088 fm for ^°Zr, and 6np = 0.201 fm 
for ^°^Pb. This type of functionals, however, has not yet been implemented 
in the RTBA model used in this study. 

In summary, the first calculation of charge-exchange spin-dipole excita- 
tions in ^^Zt and ^*^^Pb, using the particle- vibration coupling model based on 
the covariant density functional theory, is reported. Compared to the RRPA 
results, the RTBA model for particle- vibration coupling leads to pronounced 
fragmentation of the strength distribution for all three angular-momentum 
components of the spin-dipole operator. A portion of the strength is shifted 
to higher excitation energy, and the corresponding shift of the centroid energy 
is 2.5 MeV for ^°Zr, and 1 MeV for ^°^Pb. As a result of particle-vibration 
coupling a high-energy tail of the strength distribution is formed and the 
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strength in this region decreases almost hnearly with increasing energy, in 
close agreement with data. Furthermore, a model-independent SD sum rule 
has been used to determine the neutron rms radii and the thickness of the 
neutron skin. The calculated skin thickness 6np = 0.115 fm for ^"Zr and 
5np — 0-272 fm for ^°^Pb, are larger than the available empirical values, and 
reflect the high asymmetry energy of the particular relativistic energy density 
functional used in the present study. 
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